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A simple effective model of charge ordered and (or) magnetically ordered insulators is studied. The 
tight binding Hamiltonian analyzed consists of (i) the effective on-site interaction U, (ii) the intersite 
density-density interaction W and (iii) intersite magnetic exchange interaction J z (or J xy ) between 
nearest-neighbors. The intersite interaction are treated within the mean-field approximation. One 
shows that the systems considered can exhibit very interesting multicritical behaviors, including 
among others bicritical, tricritical, tetracritical and critical end points. The analysis of the model 
has been performed for an arbitrary electron concentration as well as an arbitrary chemical potential 
in the limit of strong on-site repulsion (U — > +oo). The phase diagrams obtained in such a case 
are shown to consist of at least 9 different states, including four homogenous phases: nonordered 
(NO), ferromagnetic (F), charge ordered (CO), ferrimagnetic (intermediate, I) and five types of 
phase separation: NO - NO, F - NO, F - F, CO - F, CO - I. 

PACS numbers: 71.10.Fd, 71.45.Lr, 75.30.Fv, 64.75. Gh, 71.10.Hf 



I. INTRODUCTION 

Electron charge orderings and their interplay with 
magnetism are relevant to broad range of important 
materials such as manganites, multiferroics and other 
strongly correlated electron systems [1-3]. In this report 
we discuss an effective model of charge ordered and (or) 
magnetically ordered insulators. 

We consider the following model Hamiltonian of 
a fermion lattice system in the atomic limit 

H = ^2Uh it hn - + (1) 
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where cf a denotes the creation operator of an electron 
with spin a at the site i, fii a =c~l a Ci a , hi = J^rw, 
sf = (l/2)(rij.j. — fin). U, W and J are the interaction 
parameters (on-site and intersite between the nearest 
neighbors), fi is the chemical potential, depending on 
the concentration of electrons 

n=sEi<fti>. (2) 

with n € [0,2], N is the total number of lattice sites. 
Our denotations: n Q = \{n A - n B ), n a = j| Eie a (™*)) 

m = Jf Ei(*i>. m Q = \{m A - tub), m a = §f E lGa (sf>, 
and a = A, B labels the sublattices. W = zW and J = 
z J, where z is the number of the nearest neighbors. Only 
the two-sublattice orderings on the alternate lattices are 
considered in this report. 

The interactions U,W, and J will be treated as effec- 
tive ones and be assumed to include all the possible con- 
tributions and renormalizations. One should notice that 
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in absence of external magnetic field the ferromagnetic 
( J > 0) interactions are simply mapped onto the antifer- 
romagnetic cases (J < 0) by redefining the spin direction 
on one sublattice in lattices decomposed into two inter- 
penetrating sublattices. Thus, we restrict ourselves to 
the case J > 0. 

In the analysis we have adopted a variational approach 
(VA), which treats the on-site interaction term (£7) ex- 
actly and the intersite interactions (W, J) within the 
mean-field approximation (MFA). 

For the model (1) a ground state phase diagram as a 
function of \i [4] (exact results) and special cases at T > 
(in VA) such as W = [5], J = [6-8] and n = 1 [9] have 
been investigated till now. 

Within the VA the intersite interactions are decoupled 
within the MFA, which allows us to find a free energy 
per site f(n). One can also calculate the averages: n, 
uq, m and thq, what gives a set of four self-consistent 
equations (for homogeneous phases). This set for T > 
is solved numerically and one obtains uq, m, toq and 
n (or /u) when [i (or n) is fixed. It is important to find 
a solution corresponding to the lowest energy. 

The phase separation (PS) is a state in which two do- 
mains with different electron concentration exist (coex- 
istence of two homogeneous phases). The energy of PS 
states is given by 
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f PS {n +1 n_) = kf+(n+) + (1 - fc)/_(n_), (3) 



2 



(a) 





: co . - ' 






NO 




n = 






F 




n = 1 







-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

n*J„ 



(b) 



LCO^ 

CO= 
PS2 


I=PS5 
=PS6 

F-* 


I=PS5 
=PS6 


^LCO 

CO= 
PS2 


PS4 F-» 


PS4 



0.0 



0.5 



1.0 

n 



1.5 



2.0 




FIG. 2. Finite temperature phase diagrams for W/J = —5 
(the U — > +oo limit). Dotted, solid and dashed lines indi- 
cate first order, second order and "third order" boundaries, 
respectively. Details in text. 
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FIG. 3. Finite temperature phase diagrams for W/J = and 
U — > +oo. Denotations as in Fig. 2. 



II. RESULTS AND DISCUSSION FOR U -> +oo 



FIG. 1. Ground state phase diagrams in the U — > +oo limit: 
(a) as a function of JX/ Jo and (b) as a function of n. Dotted 
lines denote discontinuous transitions (on the panel (a)). 



where E±(n±) are values of a free energy at n± corre- 
sponding to the lowest energy homogeneous solutions and 
k = n ~ is a fraction of the system with density n + . 

In the report we have used the following conven- 
tion. A second (first) order transition is a transition 
between homogeneous phases with a (dis-)continuous 
change of the order parameter at the transition temper- 
ature. A transition between homogeneous phase and PS 
state is symbolically named as a "third order" transi- 
tion. During this transition a size of one domain in the 
PS state decreases continuously to zero at the transition 
temperature. We have also distinguished a second order 
transition between two PS states, at which a continuous 
change of the order parameter in both domains takes 
place. 

In this report we present the VA results for the model 
(1) in the limit U —> +oo. This case corresponds to the 
subspace excluding the double occupancy of sites (by 
electrons for n < 1 or holes for n > 1). The obtained 
phase diagrams (for fixed n) are symmetric with respect 
to half- filling (n — 1) because of the particle-hole sym- 
metry of the Hamiltonian (1). The phases and states, 
which can exist on the phase diagrams of model consid- 
ered, are collected in Table I. 



In this section we discuss the behavior of the system in 
the limit of infinite on-site interaction (U — » +00). The 
ground state diagrams obtained in this limit are shown 
in Fig. 1. The energies of a few states are degenerated 
if the system is considered for fixed n. However, at fi- 
nite temperatures the degeneration is removed. LCO is 
the particular CO phase with uq = 0.5 and n — 0.5 or 
n = 1.5. 

The selected finite temperature phase diagrams, ob- 
tained for fixed /1 as well as for fixed n, are presented in 
Figs. 2-5. One should notice that the first order tran- 
sitions on the diagrams for fixed fi corresponds to the 
"third order" transitions and occurrence of PS sates on 
the diagrams for fixed n. The second order transition 
F-NO occurs at k B T = J \l — n\ independently of W/J. 

The labels of the critical points for phase separation are 
given corresponding to those in [10], where the Blumc- 
Emery-Griffiths model was considered. Critical points 
indicated on the phase diagrams (descriptions in corre- 
spondence to diagrams as a function of fi/ Jo, Figs. 2-5): 

(i) T - tricritical point: change of the transition order, 

(ii) C, C - isolated critical (bicritical end) points: the 
end of 1st order transition line, (iii) E, E' , E" - critical 
end points: connection of two 1st order lines and one 2nd 
order line, (iv) Z - zero-temperature critical point: the 
end of 2nd order transition line at T — 0, (v) D - zero- 
temperature highly degenerate point: the end of both 1st 
and 2nd order transitions lines at T = 0, (vi) B - bicriti- 
cal point: connection of one 1st order and two 2nd order 
lines. 
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FIG. 4. Finite temperature phase diagrams for W/ J — 0.5 
and U — > +00. Denotations as in Fig. 2. 
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FIG. 5. Finite temperature phase diagrams for W/J — 1.5 
and £7 — > +00. Denotations as in Fig. 2. 

For large negative values of W/J, a second order F-NO 
line terminates at -E-point on a first order line, which it- 
self ends at C-point (Fig. 2) . For values of W/ J — 
(Fig. 3), the transition lines meet T-point, which is con- 
nected with a change of F-NO transition order [5] . 

With increasing W/J the tricritical behavior changes 
into another one with the C"-point located inside the F 
phase (cf. Fig. 4, for W/J = 0.5). 

For W/J > 1 (Fig. 5) two new ordered phases ap- 
pear on the phase diagrams: charge ordered (CO) phase 
(riQ ^ 0) and intermediate (I) homogeneous phase, where 
magnetic order and charge order coexist (uq ^ 0, m =/= 
and ttiq 7^ 0). If W/J = 1.5 the F and CO phases are 
separated from the NO phase by two second order lines 
that connect at -B-point. The transitions between CO 



and phases with m / (F, I) are first order. The F and 
I phases are separated by second order line ending at 
£"-point. 



III. FINAL REMARKS 

It is seen that a large variety of phase transition phe- 
nomena are present in the simple model incorporating 
both magnetic and density degrees of freedom. Derived 
results are exact in the limit of infinite dimensions, where 
the MFA treatment of intersite interactions becomes the 
rigorous one. 

Our model in the limit U — > +00 is equivalent to the 
Blume-Emery-Griffiths model, which is the general spin 
5 = 1 Ising model with nearest-neighbor interactions and 
up-down symmetry. Thus obtained diagrams have the 
similar structure as those presented in Ref. 10. 

One notice that in the U — > — 00 limit of model (1) the 
magnetic orderings are completely suppressed and the 
CO phase (for W > 0) or the PS: NO-NO (for W < 0) 
state can occur. For W = only the NO is stable at any 
k B T > in this limit. 

Let us point out that in MFA the following equivalence 

occurs: 2JY,(i,j) ^ J ^( l . J ) ( s fsj + sj§r). In 
both cases the self-consistent equations have the same 
form and a magnetization along the z-axis becomes a 
magnetization in the xy-plane [11]. 
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